%x = data to be plotted
%sigma = standard devition
%mu = mean, median and mode
function f = gaussian_distribution(x, mu, s)
p1 = -.5 * ((x - mu)/s) .^ 2;
p2 = (s * sqrt(2*pi));
f = exp(p1) ./ p2; 